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ABSTRACT 


Flexural  deflections  of  several  plates  and  beams  under 
an  unknown  transverse,  concentrated,  time  dependent  force 
are  solved  for  various  edge  conditions.  The  consideration  of 
displacements  and  the  use  of  Hertz's  law  of  impact  at  the 
point  of  contact  lead  to  a non-linear  integral  equation  for 
the  contact  force  in  all  cases  of  transverse  impact.  Two 
methods  which  are  developed  in  the  previous  report  are  then 
used  to  treat  this  equation:  (a)  Generalized  Galerkin  Method; 

(b)  Collocation  Method.  Formulas  for  the  maximum  contact 
force  and  contact  time  are  obtained  which  apply  to  all  elastic 
impact  problems.  With  the  aid  of  various  curves  given  in 
this  report  calculation  can  be  shortened  to  a minimum.  Im- 
pacts on  circular  and  rectangular  plates  are  studied  in 
detail. 


1.  Introduction 


To  determine  the  deflection  and  stresses  in  a beam  or  a 
plate  struck  transversely  by  an  elastic  body  one  must  know 
the  local  compression  - contact  force  relationship  and  the 
flexural  deflection  under  an  unknown  concentrated  time  de- 
pendent contact  force.  A solution  for  the  first  is  to  use 
the  Hertz  law  of  impact  [2],  £3]#  and  [4],  which  is  known 
to  have  a wide  range  of  applicability  [5],  The  flexural 
deflections,  on  the  other  hand,  can  be  determined  by  using 
the  classical  theory  of  vibrations  of  beams  and  plates.  Con- 
sideration of  the  displacement  at  the  point  of  contact  then 
leads  to  a non-linear  integral  equation  for  the  contact  force. 
Various  authors  studied  the  problem  of  central  impact  of  a 
sphere  on  a simply  supported  beam  previously.  A discussion 
of  various  methods  given  by  these  authors  is  presented  in 
[1]. 

The  impact  problems  involving  beams  with  different 
types  of  support  conditions  and  plates  with  different  shapes 
and  non-central  impact  problems  are  as  yet  untouched,  except 
for  the  paper  by  K.  Karas  [6],  He  treats  the  problem  of 
central  impact  on  simply  supported,  rectangular  plates  using 
the  method  given  by  Lennertz  and  the  step-by-step  integration 
method  used  by  Timoshenko  for  the  beam  problem.  A criticism 
of  these  methods  may  be  found  in  [1]  also.  Briefly,  the 
method  introduced  by  Lennertz  oversimplifies  the  problem 
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leading  to  large  discrepancies  between  the  correct  and  approxi- 
mate contact  forces,  while  the  step-by-step  integration  method 
presents  a tedious  process  and  must  be  carried  out  for  each 
individual  problem,, 

In  the  present  report,  the  problem  of  non-central  impact 
on  beams  and  plates  having  general  edge  conditions  is 
formulated  in  such  a way  as  to  unify  the  transverse  impact 
problems.  Consequently,  for  the  free  transverse  oscillation 
problem  leading  to  normal  modes,  the  integral  equation  of 
impact  can  be  written  very  shortly.  Explicit  forms  for  the 
integral  equations  are  given  for  the  following  cases: 

(a)  Beams: 

1)  Simply  supported  beams 

2)  Beams  having  both  ends  clamped 

3)  Cantilever  beams 

(b)  Circular  Plates: 

1)  Simply  supported  circular  plates 

2)  Circular  plates  with  clamped  outer  edge. 

(c)  Rectangular  Plates: 

1)  Simply  supported  at  all  edges 

2)  Simply  supported  at  two  parallel  edges  and 
clamped  at  the  others. 

Two  general  methods  are  then  used  to  obtain  approximate 
solutions  of  the  integral  equations  of  contact  force: 
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(a)  Generalized  Galerkin  Method 

(b)  Collocation  Method. 

Discussion  of  these  methods  is  given  in  [1]. 

Explicit  formulas  are  given  for  the  maximum  contact  force 
and  for  the  contact  time  for  various  different  types  of  shape 
functions.  Plots  are  made  for  various  functions  to  facilitate 
computation  for  a given  problem.  In  particular,  beams, 
circular  plates  and  rectangular  plates  with  simply  supported 
edges  are  studied  in  detail. 


2.  Formulation  of  the  Problem 


A beam  or  a plate  is  struck  transversely  by  a mass  m 
having  a spherical  surface  at  the  point  of  contact,  and 
striking  velocity  vQ  , (Fig.  1) . 

The  problem  is  to  determine;  the  contact  force  F(t)  , 
deflection  w(x„y#t)  and  the  flexural  stresses. 

Deflection  and  the  stresses  are  functionals  of  the 
contact  force  F(t)  . Consequently,  F(t)  must  be  determined 
first. 

The  formulation  of  this  problem  can  be  effected  only 
under  certain  assumptions,  namely;  (a)  all  assumptions  of 
the  classical  theory  of  plates  or  beams  are  applicable; 

(b)  the  Hertz  law  of  impact  is  valid.  The  last  assumption 
states  that: 

a » k F2/3(t)  (1) 


where  a is  the  relative  approach  of  striking  bodies  and 
k is  the  Hertz  constant  [1],  [2],  [3],  The  relative 
approach  is  the  difference  between  displacements  of  the  plate 
and  the  striking  body  measured  from  the  instant  of  initial 
contact  (Fig.  1).  Hence; 


a ■ 


w 

s 


(2) 
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where  wg  is  the  displacement  of  the  sphere  under  the  action 

of  the  force  F(t)  and  wQ  is  the  deflection  of  the  plate 

at  the  point  of  contact , Here  vr  is  given  by: 

s 


wg  * vQt  - - f (t-x)  F(r)dT 
o 


(3) 


The  deflection  of  the  beam  or  the  plate  is  in  turn  ob- 
tained by  solving  the  following  differential  equation: 


DA2w  + £ « q(x,y,t) 

g dtz 


(4) 


Eh3/12  (1  - v2) 


dx‘ 


<3y 


(5) 


Here  E is  Young's  modulus,  h is  the  thickness,  p is 
the  weight  density  per  unit  area,  g is  the  acceleration 
due  to  gravity,  q is  the  load  per  unit  area,  Equation  (4) 
reduces  to  beam  equation  if  one  takes  El  for  D and 
y ■ d/dy  ■ 0.  In  this  case,  p and  q must  be  interpreted 
as  the  weight  density  and  the  load  per  unit  length,,  In  the 
following  analysis  beam  and  plate  problems  are  interpreted 
as  one.  The  difference  is  in  the  additional  dimension  y 
and  the  interpretation  of  D and  p which  can  be  adjusted 
easily  to  the  method  of  solution.  The  resulting  formulas, 
however,  will  be  given  separately  for  beams  and  plates. 

The  concentrated  load  F ( t ) is  obtained  as  a limiting 
process  of  a uniformly  distributed  load  q over  a small 
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o 

area  e having  the  shape  of  the  plate  boundary.  In  beams, 
q is  distributed  uniformly  over  a range  e . Thus,  by 
letting  q— * oo  while  e— >0,  the  contact  force  F(t)  is 
obtained: 


F (t) 


lim  f 

e — * 0 


q <3A 
(A) 


(6) 


q -►  oo 

Hence  equation  (4)  with  q subject  to  (6)  must  be  solved 
under  initial  and  boundary  conditions  to  obtain  wQ  . Once 
this  is  done,  the  force  F(t)  can  then  be  evaluated  by 
solving  the  integral  equation  obtained  by  combining  equations 
(1)  to  (3),  namely: 


kF2/3(t) 


v t - — f 

O m o 


(t“T)  F(t)  dT  * W 


(7) 


i 


! 


I 

I 


i 


} 


V 


3.  Flexural  Deflection 


Deflection  w of  a plate  or  a beam  is  uniquely  deter- 
mined when  equation  (4)  is  solved  under  appropriate  initial 
and  boundary  conditions.  Initial  conditions  are: 

w(x,y,0)  - 0 w(x,y,0)  - 0 (8) 

where  the  dot  represents  differentiation  with  respect  to 
time.  Boundary  conditions  used  in  practice  have  great 
variety.  These  conditions  consist  of  specifying  any  two 
of  the  following  four  quantities  or  their  linear  combinations 
on  the  boundary: 

Given:  w , dw/dn  - Ifl  . V (9) 

n 9 n 

Here  n represents  the  external  normal • to  the  boundary  curve 
and  is  the  bending  moment  and  V is  the  transverse 

shear  resultant  which  is  represented  by  the  Kirchhoff  condition 
in  the  classical  plate  theory. 

In  what  follows,  first  a general  solution  of  equation  (4) 
is  found  which  satisfies  initial  conditions  (8)  and  some 
arbitrary  boundary  conditions, 

Dirac-delta  function  8 (x  - x ) is  defined  by: 

o 

(oo  for  x ■ x 

and  / 5 (x-x  )8  (y-y  )dA  - 1 

0 for  x f xQ  (A)  u ; 

(10) 
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where  the  double  integral  in  equation  (10)  is  taken  over  the 
whole  area  A of  the  plate.  After  Schwartz's  book  [7], 
the  old  argument  about  the  existence  of  such  a function  can 
be  looked  upon  as  a solved  problem. 

Let  the  time  dependent  concentrated  load  F(t)  be 
applied  to  the  point  with  coordinate  (xQ,  y ) . Differential 
equation  (4)  can  be  written  in  a convenient  form  by  using 
equations  (6)  and  (10) : 

DA2w  + (p/g)  w - 8 (x-xQ)  8 (y-yQ)  F(t)  (11) 

The  solution  of  this  equation  is  a linear,  homogeneous 
functional  of  F(t).  Thus: 

t 

w(x,y, t)  ■ f G(x,y, t-x)  F(x)  dx  (12) 

o 

Substitution  of  equation  (12)  into  equations  (11)  and  (8) 
gives : 

DA2G  + (p/g)  0 - 0 (13) 


G (x,y,0)  - 0 , (p/g) G(x,y,0)  * S(x-xQ)  S(y-yQ)  (14) 


The  general  solution  of  equation  (13)  satisfying  the  first 
initial  condition  given  by  the  first  of  equations (14)  and  all 
boundary  conditions  is: 


m,n 


a g 
mn  *mn 


U,y) 


sin  a t 

uni 


(m,n  * 1| 2| • • • ) 


k 


G 


(15) 
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where  a are  arbitrary,  and  g„„(x.y)  are  free  of  arbi- 

trary  constants  and  satisfy  all  boundary  conditions  and  the 
differential  equation; 

A2gffln  " Pmn4gmn  - 0 , (16)  p 4 - Pa_n^/gD  (17) 

mn  mn  mn  ° ran  mn 

Functions  gmn  are  the  eigenfunctions.  Quantities  are 

the  eigenvalues  and  are  determined  from  a frequency  equation 

which  will  be  known  when  the  boundary  conditions  are  explicitly 

given.  a „ must  now  be  determined  from  the  second  of 
mn 

equations  (14) ; 

s (p/g)amnamngmn(x„y)  - 6(X-X0)  8(y-y0)  (18) 

m,n 

When  the  eigenfunctions  are  not  orthogonal  this  represents  a 
difficult  problem  and  there  is  no  universal  method  of  deter- 
mining the  amn  o In  many  cases  it  is  possible  to  use  the 
method  used  for  orthogonal  eigenfunctions  except  that  the 
series  (15)  must  be  employed  to  sum  the  resulting  series. 

In  many  important  practical  problems  the  eigenfunctions  are 
orthogonal.  Consequently,  affin  can  be  determined  simply  by 
multiplying  equation  (18)  by  g„e  , integrating  over  the 
total  area  of  the  plate  and  considering  the  condition  of 
orthogonality; 
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S gm„g, 

(A)  mn  1 


Hence: 


amn  * 'WWo1  «'•>  bmna., 


Consequently: 


w(x,y,t)  = (g/p)  mEn  (l/bmn)gmnU0,y0) 


/ (1/a  ) F (t) sin  a_n(t-x)dx 

_ mn  mn 

o 


The  deflection  at  the  point  of  contact  is  therefore  found  to 


w » (2/M)  £ c S (l/a__)  F(t)  sin  a (t-x)  dx  (22) 
o _ _ mn  “ mn  mn 


m,n  o 


°mn  ' 9nm  ^o^o1  A'2bmn 


where  A is  the  total  area  in  case  of  plates  and  total  length 

in  case  of  beams,  and  M is  the  total  mass. 

When  the  free  vibration  problems  of  beams  or  plates  are 

solved  the  eigenfunctions  g and  eigenvalues  p 

mn  mn 

(consequently  ) will  be  known.  Thus,  the  determination 

of  b__  from  equation  (19)  and  <2  from  (23)  presents  no 
mn  mn 

difficulty. 
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A few  practical  examples  for  beams  and  plates  are  given 
below.  In  all  examples  given  below,  the  eigenfunctions  are 
represented  by  gn  or  g^  and  the  frequency  equation  is 
marked  by  (F.E.).  Symbols  n and  m represent  positive 
integers. 

a)  Beams.  Eigenfunctions  and  eigenvalues  of  the  follow- 
ing examples  are  known  [8],  Thus,  bR  and  cn  are  cal- 
culated by  using  Equations  (19)  and  (23)  above,  since  in  all 
cases  the  eigenfunctions  are  orthogonal.  The  origin  is  at 
one  end  of  the  beam. 

1)  Simply  supported  beams 

gR(x)  * sin  nnx/L  # an  “ n^n^g  EI/p  (F.E.)  (24) 

2 

b - L/2  , c * sin  nnx/L 

n ' ' n o' 


gn(*) 


2)  Beams  having  both  ends  clamped 
sinh£pn (2x-L) 


sinhipnL 


sin£pn(2x-L) 

sinjpnL 


cos  pnL  ■ sech  pnL  (F.E.) 

1 rPnL  ♦ sinh  pnL  pnL  + »in  P„L 
bn  pn  t i-oosh  pnL 


t , 1 

1-cos  p L J 


} (25) 


c - g (x  )L/2  b 
n ^n  o n 
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3)  Cantilever  beams  clamped  at  x ■ 0 

cosh  p x - cos  p x sinh  p x - sin  p x ^ 

a (x\  . n n „ n *n 

ynv  cosh  pnL  + cos  pRL  sinh  pnL  + sin  p^L  • 

cos  p L - - sech  p L (F.E.) 
n n 

-I  -2 

bn  * pn  {cosh  pnL  ” cos  pnL)  (pnL 

c » g 2 (x  ) L/2b 
n 3n  o n 

4)  Beams  pinned  at  x ■ 0 and  clamped  at  x. 

sinh  p x sin  p x 

^n  sinh  pnL  sin  p^L  ' 

tan  pnL  - tanh  pnL  (F.E.) 

t>n  * 7 (sin  2 PnL  “ sinh”2  pRL)  , 

°n  * «n2(xo>  L'2\ 

b)  Circular  Plates.  Eigenfunctions  and  eigenvalues 
of  the  following  examples  are  known  [8].  Equations  (19)  and 
(23)  above  give  bn  and  cn,  respectively,  since  in  both  of 
these  cases  the  eigenfunctions  are  orthogonal.  The  origin  is 
taken  at  the  center  of  the  plate.  The  outer  radius  is  a . 
Polar  coordinates  r and  0 are  used. 


>(27) 


^ (26) 


cos  p L - cosh  p L 
n n 

sinh  p L 
n 


sin  p L + 
n 


J 
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1)  Simply  supported  circular  plate 
9mn(r>e)  - tJn(pmr)  + B Vpiiir)} sln  (n6  + V ' 

(28) 

B“  * Jnlpma)/In(pma> 

where  J and  I are  the  Bessel  functions  of  the  first 
n n 

kind  with  real  and  imaginary  arguments. 


lpma)  „ Jn»l  lpmai 
an  tpmat  Jn  lpma) 


2pma 

X-v*" 


(F.E. ) 


(29) 


Here  v is  the 

b * 2na2 
mn 


c “ g 
mn  ymn 


Poisson's  ratio. 

j 2(p  a)  r1  * v Vp  - 

n *n  L v - 1 

(r  , na2/2b 
o'  o mn 


(30) 


2)  Circular  plate  clamped  at  the  outer  edge 

g (r,©)  - [Jn<Pmr)  + B In(Pmr^sin  + 

3nin  n m n m o 


B-  ‘ Jn(pma)/In(pma) 


Jn+1  (pma) 


+ 


Vl  (pma). 


(31) 


0 


(F.E.) 


(32) 
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mn 


S™2  Jr.2  <Pma) 


mn 


g__2  /4Jn2  (Pma) 

’mn  o o n m 


(33) 


c)  Rectangular  Plates.  For  rectangular  plates  a quick 
method  similar  to  the  one  used  for  beams  can  be  used  to 
obtain  boundary  conditions  which  will  produce  orthogonal 
eigenfunctions. 

The  following  is  constructed  from  Equation  (16) : 


b a 


b a 


/ r g A g - g A g )dxdy  * (p  - p )f  / g g dxdy 

j j 'yrs  ymn  ymn  yrs'  y v*-mn  ^rs  u j ymnyrs  y 
o o o o 

(34) 


Integration  by  parts  of  the  left  side  of  this  equation  leads 
to: 

b a 


, _3  , ,2 

a a v a DC  g og^„  o g 

(p  4-p  4)  // g g dxdy  - / [g  — 25- Sr1 

rs  o o mn  rs  o rs  djt5  ax2 


3 grs  a9mn 

"S*"35" 


.3 
o g 

us  . * 

axs  mn  17 


(35) 

X ■ 

dy*dx  x * 0 


a 3gmn  3 ^ * * ady 


and  a similar  equation  for  y * 0 and  b , which  may  be 
obtained  by  interchanging  x with  y . The  bracket  on  the 
right  side  of  equation  (35)  is  zero  for  two  types  of  practical 
boundaries  --  namely,  simple  support  and  clamped  edge.  Con- 
sequently, these  two  types  of  support  conditions  in  any  order 
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lead  to  normal  modes.  Many  other  conditions  can  be  obtained 
by  making  the  bracket  on  the  right  side  of  Equation  (35) 
zero,  with  the  proper  choice  of  g and  its  derivatives. 

The  following  two  examples  are  most  practical: 

1)  Rectangular  plates  simply  supported  at  all  edges. 
Two  perpendicular  sides  with  lengths  a and  b are  taken  as 
the  coordinate  axes,  origin  being  at  a corner: 


/ \ nmx  nny  2 gD  4 

g__  (x,y)  * sin  — — sin  < * , a m P__ 
mn  a p ' mn  p mn 


pmn2  " H**/*)2  * inn  lb) 2]  , b^  - ab/4. 


mn 


nnyrt  9 

2 (sin  — sin  - i— -°-)  2 


mnx. 
"a" 


(F.E. ) 

i(36) 


2 ) Rectangular  plate  simply  supported  along  x « 0.  a ; 
clamped  along  y ■ 0.  b . 

^mn^'Y^  “ cos^  ay  " cosPy  - B(sinPy  - sinh  ay) 

B ■ (cosh  ab  - cosPb)  / (cosfib  - ^ sinh  ab)  >(37) 

a * [p2  ♦ (nm/a)2]^  , P * [p2  - (mn/a)2]^ 

(cosPb  - cosh  ab)2  + (sinPb  - sinh  ab) 

a 

(sinpb  ♦ sinh  ab)  ■ 0 (F.E.) 


(38) 
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fa  » iab(l  ♦ (sinPb  - sinh  ab)"1 

mn  a a 

a 

2 

• f [ (1  - ) sinh  ab  sin  fib  - 2 ^ cosh  ab  cos  fib] 

<■  a 

(cosh  ab  + cos  ab)  ♦ 2 ^ (cosh  ab  + cos  fib)}  , 


c * g__(x,y)  ab/ 2b 
mn  mn  o o mn 


4.  The  Integral  Equation 


In  the  previous  chapter  it  is  shown  that  in  all  impact 
problems  the  flexural  deflection  at  the  point  of  contact  is 
given  by  Equation  (22).  The  integral  equation  (7)  of  impact 
can  be  written  in  a common  non-dimensional  form  as: 


L(f ) * af2/3(x)  - x + / K(x-y)-  f (y)dy  - 0 


(40) 


f(x)  « T0F(TQx)/mv^  , K(x)  - x + 2n  E JBL  sinwm„x  , 


c 

_i 

ea- 


rn, n mn 


mn 


H - m/M 


a)  * a T 
mn  mn  o 


t/T 


k m2'\  -1'3  T'5'3 


> (41) 


where  T may  be  chosen  in  a suitable  manner.  There  c „ 
o mn 

and  m.  are  given  by  the  solution  of  the  flexure  problem, 
mn 

Another  form  of  Equation  (40) , which  will  be  used 
later,  is: 


a0(x)  - x ♦ S K(x-y)  03^2(y)  dy  - 0 
o 

$)(x)  - f2/3(x) 


(42) 
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The  Generalized  Galerkin  method  is  developed  in  [1]. 

This  method  gives  an  approximate  solution  for  any  problem 
having  an  equation: 

L(f ) - 0 (43) 

Let  a function  f which  depends  on  the  arbitrary  independent 
infinitely  many  parameters  (a^#^, . ..)  and  independent 
variable  x,  be  chosen  such  that  either  df/<3a^  satisfy  the 
boundary  conditions  at  x « x^  and  x = x2  , imposed  on 
equation  (43)  or  L(f)  = 0 there.  This  function  f(x„  apa^,...) 
obviously  will  not  satisfy  equation  (43) . This  function 
represents  a solution  of  equation  (43)  if  a^  are  solutions  of 
the  following  set  of  non-linear  algebraic  equations: 


x2 

S L[f  (xta1<ra2,  „ . „ ) ] dx  - 0 
xi  1 


(44) 


When  i is  finite  f (x,a^,a2, . . , ,a^)  is  an  approximate 
solution. 

Various  functions  of  the  following  type  are  suitable  to 
use  as  f (x,a1,a2, o . . ) : 

(a)  f « a^  sin  a2  x 

2 

(b)  f * a^sin  a2x 
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(c)  f • a^  sin^^  agX 


N-l 

(d)  f • E a.,  sin  i aMx 
i=l  1 N 


Hence  equation  (44)  takes  the  form: 

/n/&N  L(f ) (df/3ai)dx  - 0 (i  - 1,2,... ,N)  (45) 


Note  that  df/da^  and  df/dag  for  cases  (b)  and  (c)  satisfy 
the  end  conditions  at  x - 0 , n/ ag  . However,  in  case  (a), 
one  end  condition  is  violated  since  df/dag  does  not  vanish 
at  x • nfa-2  , end  of  the  contact.  Nevertheless,  in  all 
examples  worked  out,  very  good  approximations  are  obtained. 
Cases  (b)  and  (c)  give  better  results  than  case  (a).  (See 
Figures  18,19  , and  20.)  In  case  (d)  ajj  is  chosen  the  same 
as  a2  of  case  (a)  since  ag  of  case  (a)  is  found  to  be 
nearly  equal  to  its  exact  value. 

Explicit  forms  of  Equation  (45)  are  given  below  for 
cases  (a),  (b) , and  (d) . Computations  are  carried  out  for 
cases  (a),  (c) , and  (d). 


Case  (a) 

(£  + 2\i  a^)  4 d.y  • * 1 

<£-£nS)  *!  + ■&<!  -y^2'3  - 1 


> (46) 
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Elimination  of  X^^  between  these  two  equations  leads  to: 


Ax  - 28 /f 29  - n (10B  + 24^)] 


y * a a «^3/rc 


15  - 10 ii  (B  + 8^) 
d 282/3[29  - n (10S  + 24ax)] 


(47) 


Case  (b) 

<F-  W + *1+  d'^l2/3  * 1 


2/3 


(48) 


Similarly,  X-^3  can  be  eliminated  from  these  equations 
to  obtain  equations  resembling  equations (47) . 

Case  (d) 

E (1  ♦ 2n0,.)  X ♦ [1  - 4 (-l)j  + 1]  X.  ♦ ^ia,  • X. 

i.ifj  1J  i ^ J J J 


+ (-l)j+1  YH j X^'"  - 1 (j  - 1,2,...) 


2/3 


(49) 
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where : 


f>ii  + t'1’1  7 


S c A (Tv*  ) f 
_ mn  mn 
m#n 


8 ' h * «1  ' £.  °mn  DSm> 


y - a a25/3/n  , d - n1/2  T (4/3)//“ {11/3)  * 1.68257 


81  ‘ z 0onDl  (\»n>  « *1  ' 2 °mnBi  Sm1  > 

m,n  m,n 


^ ‘ m!n  0nnBij  <5W’ 


B1  ^mn1  ’ (1  •>’mn2)‘2  cos^mn  + (1  ‘ \in2)‘3 


* «1  +0(*\»,‘1  Sinrt>-mn 

Bi  Kn'  - t1  ' ^/i»23'1 


Bi;j<Snn>  -t(1  * <Snn'i>2]'1  C1'  KJ^'1 


'"Nan1"  sinn\m 


* VaN  * V-«  ““n,n/aM  " 1,2,...) 


mn  mn'  N 


f 
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“1  ‘ A °“>A  ‘V1  ' 8 " A °nmD  (\m> 

m,n  m,n 

a'  (^n)  - 2(2 /Xmn) 2 - (1  - ^V1  - i (^.mn)*3 
‘ <!  - KA'2  slnnv  . 


* <2'Sm>2  + T'1  - K2''1  'TO'3 

2 1 2-2  l (51) 

• sin n\-  (2/\mn)  (1  * 4*.mn  ) * cosn\ 

mn  mn  4 nin  mn 

d'  - (14/5ti1/2)  r(7/6)  / r(5/3) 


K'-,2^  3t]  . «/  (Z  Isin  sin  judu  (i,j  * 1,2,...(N) 

J o i 

(52) 


2/3 


For  N * 3,  have  the  explicit  forms  given  below  which  are 

obtained  by  using  the  binomial  expansion. 

t|1  - 1.682621  + 0.101974  (3  Eg /E^  > 

t)2  - 1.223724  (2  E^)  V (53) 

rjg  - 0.152967  ♦ 1.265714  (3  Eg/E^ 
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A<  (\  ) • D(\)  » B,  (V)  f B,  . (V  ) , (i  f j)  f 

l mn  ' mn  i mn  lj  mn 

» t 

A (^mn)  * and  D (^mn)  are  plotted  on  Figures  2 to  5.  Hence 
for  any  impact  problem  these  curves  can  be  used  to  sum  up 
the  series  involved  in  the  expressions  of  and  5 . 

If  the  contact  time  nf a^  is  chosen  by  an  approximate  method 
as  suggested  by  Zener  and  Feshbach  and  E.  H.  Lee  [1]  or 
simply  taken  as  the  Hertz  time  (as  this  seems  to  give  a better 
approximation)  then  these  series  can  be  computed.  One  can 
then,  for  example  in  case  (a),  use  the  first  of  equation  (46) 
to  determine  the  maximum  contact  force  7^  . This,  of  course, 

assumes  that  the  second  equation  of  (46)  is  satisfied.  Since 
this  will  not  be  the  case  in  practice,  this  reference  contact 
time  rt/a2  may  be  used  on  the  right  sides  of  equation  (46) 
to  obtain  a better  a2  and  • 

For  particular  beams  and  plates  this  process  may  be 
further  shortened  by  giving  curves  for  » ai  » 

In  this  report  this  is  done  for  the  following  simply  supported 
plate  problems  under  central  impact: 

(i)  Circular  plate.  Fig.  6 

(ii)  Square  plate.  Fig.  7 

(ill)  Rectangular  plate,  with  a/b  * 2,  Fig.  7 
(iv)  Rectangular  plate  with  a/b  ■ 4,  Fig.  7 


Similar  curves  were  given  in  the  previous  report  for  simply 
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sup  ported  beams  under  central  impact  (Fig.  20,  [1]).  For 
other  types  of  beam  and  plate  problems  this  may  be  done  once 
and  for  all,  reducing  the  problem  of  impact  to  a routine 
one.  Further  discussion  and  application  of  this  method  to 
various  problems  are  given  in  Chapter  8. 


6 . Collocation  Method 


This  method  consists  of  satisfying  the  integral  equation 
(40)  or  (42)  at  various  characteristic  points.  In  a contact 
force  curve  the  important  points  are:  origin,  the  point  of 

maximum,  the  end  of  the  contact  time.  The  contact  force 
curve  differs  slightly  from  a symmetric  shape.  Thus,  choice 
of  a symmetric  curve  is  reasonable.  Let  x be  the  un- 

C 

known  non-dimensional  contact  time.  Also  let 


f (x)  - A3X(x)  , X(0)  - X(x  ) - 0 , X(x  /2)  - 1 

G G 

X (54) 

S(x)  - f K(x-y)  X(y)dy 
o 


3 

where  A is  the  unknown  maximum  aplitude.  Equation  (40) 

must  be  satisfied  at  x - xq/2  . Also  df/dx  obtained  from 

Equation  (40)  must  be  zero  at  this  point.  This  leads  to  two 

equations  to  evaluate  A and  x , namely: 

c 


A3  - S_x_1(xc/2) 

V2  " a S,x*2/S  (xc/2)  * S(x0/2)  S x-1(x0/2) 


(55) 


where: 


26 
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Equations  (55)  are  also  valid  if  one  uses  Equation  (42)  and 
choses: 


0(x)  - A2Y(x)  , Y (0)  - Y(x  ) - 0 # Y (x  /2)  - 1 

c c 


S(x)  ■=  /X  K(x-y)  Y3/2(y)  dy 
o 


(56) 


Hence,  if  a shape  function  x(x)  or  Y(x)  , having  the 

conditions  stated  in  Equation  (54)  or  (56),  is  chosen,  A 

and  x can  be  calculated  from  Equation  (55). 
c 

Three  practical  shape  functions  are: 

o 

(a)  X(x)  • sin  tix/x  , (b)  X(x)  **  sin  nx/x  , 

C C 

(c)  Y(x)  * sin  nx/x„ 

c 

The  results  for  S and  S in  these  cases  may  be  written: 

# X 

5(xc/2)  » (xc/7i)2(d1  + ^iR1)  , 

S,x(xc/2)  = (xc/n) (d2  + 2uR2)  (57) 


where: 
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Case  (a) 


dl  “ Y " 1 ' d2  * 1 ' Ri  “ mZn  cmnpi(>y*mn)  (i  “ 1,2) 


2v-l  * -1 


" (1  ^m‘)_A  sin  ‘ 1} 

p2(\in)  “ (1  "x'mn2)  1 cos  2*%n 


> (58) 


^mn  "“in  xc^n 


Case  (b) 


* (n  - 4)/16  , dg  “ n/4 


1 


Pl<\n„>  “ £4V,'2  - (1  * Kn2’’1  (1  * 4\nn*2  008  T’Sm’1 


p2(\m>  -^mn*1(1  * Kn^'1  sin 


V(59) 


Case  (c) 

dx  - 0.457906  , Rx  - 0.915310  - O.OIISOORjl 

d2  - 0.811409  , R2  * 0.915310  Rg#  - 0.033901R2 

i > . . 

R^  * same  as  R^  of  case  (a)  \ Rg  * same  as  R2  of  case  (a) 

R1  ” m?n  °mnpl  ^mn1'  R2  “ m^n  °mnp2  ixmn)  >(60) 

’ Cl  - C(3^n)  sin  Xn  + « 


PE'^mn1  ‘Cl  * fc-an/3’2)'1  008  Hm 


J 
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Quantities  P^^)  # P2^mn*'  P1  ' n2*  v\m^  9iven  by  equations 

(58)  to  (60)  are  plotted  as  functions  of  \ in  Figures  8 to  12 . 

mn 

Consequently,  for  any  type  of  plate  problems,  series  R^  and  Rg 
can  be  obtained  by  summing  the  individual  terms  given  by  these 
figures.  In  particular,  these  series  summed  for  simply- 
suooorted  circular  elates  and  simply-supported  rectangular  plates 
with  the  side  ratios  a/b  =1,  2 and  4.  Consequently,  for  the 
latter  plates  under  central  impact  we  can  read  the  values  of  R^, 

Rg,  and  Rg'  from  Figures  13  to  17,  or  from  Tables  4,  5 and 
6.  Hence,  in  commutation,  we  first  select  a contact  time  x 

c 

then  compute  from  the  last  of  equation  (58).  Next  we  read 

R-^  and  R2  from  these  figures.  Using  equation  (57)  we 

calculate  the  right  side  of  the  second  of  equations  (55).  If 

this  comes  out  to  be  x / 2 which  we  chose  to  start  with,  the 

c 

3 

calculation  is  finished.  A is  then  calculated  by  the  first  of 

equation  (55).  If  it  does  not  we  try  another  x / 2 value. 

c 

Usually,  two  or  three  trials  are  sufficient  to  determine  x 

v 

satisfying  the  second  of  Equations  (55).  This  method  is  faster 

than  the  generalized  Galerkin  method.  It  is  not,  however,  as 

general  and  consequently,  in  some  extreme  cases  of  impact  problems, 

such  as  multiple  imoact,  it  may  not  be  as  satisfactory  as  the 

other  method.  Comoutations  are  made  for  circular  and  rectangular 

plates  using  forcing  function  shaoes  given  by  cases  (a)  and  (c). 

The  results  are  dotted  on  figures  18  to  20  next  to  those  obtained 

■% 

by  the  Generalized  Galerkin  method. 
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In  the  case  of  square  olates,  the  result  of  K.  Karas  [6] 
obtained  by  the  reliable  stepwise  integration  method  is  plotted 
in  Fig.  19.  An  examination  of  these  curves  indicate  that  the 
approximat ion  offered  to  the  impact  problems  is  satisfactory 
for  all  engineering  purposes.  Further  examole  calculations  are 
given  in  Chapter  8. 
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7.  Deflection  and  Flexural  Stresses 
Deflection  w can  be  calculated  from  equation  (21)  as  soon 
as  F(t)  is  determined  by  one  of  the  previous  methods. 

All  equations  obtained  for  contact  force  can  be  put  into  the 

form: 


F (t)  = Z 

i = 1,2 

F (t ) = 0 


F^sin\^t 


0 £ t £ T 

c 

t i T 


(61) 


where  F^  and  are  known  or  can  be  determined  from  the 
expressions  of  the  contact  force. 

Substitution  of  the  Equation  (61)  into  equation  (21) 
leads  to: 

w(x,y,t)  * (2/M)Z  cmn  [g(x,y)/g(xo,yo)]  (F± / X^) 
m,n,i 

[1  - (a  . /X. ) 2]  1 [ (X. /a  ) sina  t - sinX.t], 

^ mn  i u l mn  mn  l 


For  0 4 t * T 


(62) 


w(x, y , t ) = (2/M)Z  c Jg(x,y)/g(x  y )][Bm  sina  t 
• mnL ’ ^ o' J o JL  mn  mn 

m,  n 


♦ 2 cosat] , For  t i T 
mn  mn  c 


(63) 


where : 


3 « Z (F ./  [1  - (a  /X. ) 21_1  [l-(-l)1  cosa  T ] (X.  /a  ) 

mn  ^ l l u mn  l 1 L mn  c J l mn 


mn 


Z (F./  X.2)[l  - (a  /Xj  ) 21_1  (-1) 1 (X.  /a  ) sina  T 
^ l l u mn  i 1 l mn  mn  c 


) 


(64) 


The  deflection  at  the  point  of  contact  is  obtained  by  taking 


x = xQ,  y = y in  Equations  (62)  and  (63) . 
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The  flexural  stresses  a , a and  shear  stresses  t 

x y xz 

and  t can  be  determined  from: 
yz 


°x  - «Vh‘ 


- ’xz  * 3V2h 


(65) 


°y  * 6My/h‘ 


• Tya  - 3Qy/2h 


where ; 


M = -D (w  + vw  ) , M = -D(w  + vw  ) 
y ,xx  , yy  y ,yy  ,xx 


Mxy  ■ 'Myx  = ‘D(1-',)w,xy 


«x  = -Dlw.xxx  + w,yyx>  • Qy  * 'D  (w.xxy  + w.yyx> 


) (66) 


Here  M , M are  bending  moments,  M twisting  moment, 
x y xy 

Q and  Q are  the  transverse  shear  forces,  all  per  ur.it  length, 
x y 

The  flexural  stresses  contain  second  and  third  derivatives  of 

deflection,  consequently,  convergence  is  very  slow.  As  a matter 

of  fact,  x and  t diverge  at  the  point  of  contact, 
xy  yz 

This  fact  is  well  known  in  plate  theory.  Bending  stresses 
a and  a are  slowly  convergent  everywhere.  This  is,  of 

x y 

course,  the  case  even  for  the  static  deflection  under  a concentrated 
load. 

As  pointed  out  in  the  previous  reoort,  [1],  the  shaoe  of 
the  contact  force  has  little  influence  on  the  deflection  of  the 
plate.  As  a matter  of  fact,  if  one  considers  the  extreme  case 
of  the  Dirac  5- function  type  contact  force  having  the  same  impulse 
as  F(t)  and  olace  it  at  the  half  of  the  contact  time  interval, 
one  finds  that  the  deflection  oroducedthis  way  differs  little 
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from  the  actual  deflection.  Fig.  22  . Of  course,  a shape  closer 
to  the  actual  shape  will  give  better  results.  Therefore,  the 
maximum  error  in  deflection  would  be  reached  by  considering  a 
8-function  type  of  contact  force.  This  fact  also  indicates  that 
the  deflections  cannot  be  used  as  a criterion  of  accuracy  for  the 
shape  of  F (t) . The  situation  for  stresses  is,  of  course, 
different . 

Deflection  at  the  centers  of  a simply-succorted  square 
plate  and  of  a circular  plate  due  to  a central  impact  are 
computed.  The  result  is  clotted  in  Figures  21  and  22.  The 
result  for  a 8-function  type  contact  force  for  the  case  of  the 
square  plate  is  also  clotted  in  Fig.  22.  Deflection  given  by 
K.  Karas  [6],  based  on  the  reliable . step-by-step  integration 
method  which  is  also  drawn  in  Fig.  22,  compares  very  well  with 
these  approximate  solutions. 
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8 . Illustrative  Examples 

Calculations  for  the  case  of  central  impact  on  simply-support- 
ed beams  are  given  in  our  previous  reoort.  In  the  present 
report,  examples  are  worked  out  for  the  case  of  central  impact  on 

(i)  Simply-suo ported  circular  plates 

(ii)  Simply-supported  rectangular  olates 
having  size  ratios  a/b  ■ 1,  2 and  4. 

Both  the  Generalized  Galerkin  method  as  well  as  the 
Collocation  method  are  used  to  determine  the  contact  force 
and  contact  time.  In  all  problems  the  striking  body  is  assumed 
to  be  a sphere.  This,  of  course,  is  immaterial  as  far  as  the 
theory  is  concerned,  since  all  we  need  is  to  require  that 
the  striking  body  should  have  a known  radius  of  curvature  at 
the  point  of  contact,  in  accordance  with  the  Hertz  theory  of 
impact.  Once  this  radius  of  curvature  is  known,  the  Hertz 
coefficient  k can  be  determined  [4"].  For  a sphere  striking 
a elate  the  equation  for  k is  -particularly  simple  and  is  given 
by: 

k = 1.230520  (E2r)"1/3  (67) 

s 

where  E is  the  common  moduli  of  elasticity  of  sphere  and  plate 
rg  is  the  radius  of  the  sphere.  In  all  the  calculations, 
the  following  quantities  are  taken  the  same: 

E = 2.2  x 10^  kg /cm2  , v = 0.3 

rs  = 1.00  cm  , h = 0.8  cm 

vQ  = 100  cm/sec  , y * C. 00796  kg/cm3 

where  v is  the  common  Poissorfs  ratio  of  both  plates  and 
spheres,  h is  the  plate  thickness  and  y is  the  common 
weight  density  of  the  plate  and  the  sphere.  Hence  k,  for 
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0.727456  x 10'4  cm/kg2/3 


the  examples  of  this  report,  is  given  by: 
k 

(i)  Circular  Plate  simply-supported  at  the  outer  edge: 
Radius  of  the  plate  = a = 10cm. 

Values  of  o a which  are  the  roots  of  the  frequency 
* m 


equation  (29)  are  tabulated  in  Table  1.  am  and  c^  are 
calculated  from  equations  (17)  and  (30)  where  for  the  case 


of  central  imoact  in  expression  (28)  of  g „ (r,6)  we  must  take 

mn 

n * 0 and  y = it/ 2.  Hence  for  this  case  we  have: 
o 


- - 1 r 1 - 1 

-a  T C' X(5"aT  I (dT) 

o m o *m 

. „Pma  Jl(-°®a)  ,-l 

2T^  J„(pa)  J 


1C 


1 + v 


pa  9 

2^)  2 


> (68) 


am  = (gD/p) 1/2  pm 


J 


Values  of  c are  giver,  in  Table  1 for  v =0.3. 

In  what  follows  the  contact  force  and  the  contact  time 
are  given  in  units  of  kilograms  and  seconds. 

(a)  Typical  Calculation  for  Generalized  Galerkin  Method 
We  use  shape  (a),  that  is,  we  assume  that  f(x)  = a^sina2x 
1st  approximation:  Take  jjl  ■ 0,  then  equation  (47) 

gives: 


^ = 0.965517  , y = 0.314688 

2nd  approximation: 

V1  a2  » (n  Y/a)3^5  = 4.64944  x 1C4,  hence 


**1G  = to10^a2  “ 0. 066897  1 1 — ) (Fig.  6)  ) 
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* 5.35  , 8 = -46.8 

Hence,  equation  (47)  gives: 

7^1  * 0.837186,  To'1a2  = 4.6127  x 104,  \1Q  = C.0674 

Y = 0.310553 

Continue  in  this  way  until  the  difference  between  two  subsequent 
values  of  y or  or  both  are  small  enough.  Even  above, 
the  second  approximation  gives  a very  satisfactory  result. 

The  values  of  the  contact  force  and  time  are  given  by: 

F ( t)  = F1sin  \t  , ?±  « 131.509  kg 

\ * 4.6127  x 104sec  ^ , T = C.681  x 10  4sec. 

c 

where  Tc  is  the  contact  time.  The  result  is  plotted  on 
Fig.  18. 

(b)  Generalized  Galerkin  Method  having  shape  (d) 

F(t)  = IE  F.  sin\t,  shape  (d) . 
j-1,2.3  J 

Three  equations  given  by  equation  (49)  when  solved  give: 

Fx  - 133 ; 427  kg  , F2  = 1.506  kg  , F3  - -11.981  kg 

Fmax  = 14 5 • 408  k<?  . Tc  = 0.679  x lO-4  sec 

In  this  calculation,  \ of  the  foregoing  calculation  is  used. 

The  result  is  plotted  in  Fig.  18. 

(c)  Typical  Calculation  for  the  Collocation  Method 
We  use  shape  (c) , that  is: 

F(t)  = A3  sin  3/2nt/Tc 

Equation  (55),  upon  substitution  of  equation  (56),  gives: 
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A3  - nmv0/T0  (dj  - %tHj) 


A 


k3m2v0"17i2  (d2  + ^xR 2 ) 


(69) 


1st  approximation: 


A = 176.863  kg 


2nd  approximation: 
**10  = a10Tc^u  = C* 


Take  ^ = 0,  then  Equation  (69)  gives: 
, Tc  = 0.684964  x 10'4sec. 
Calculate 

067815  » Fig.  13  >RX  = 6.86, 


R2  * 8.42 

Substitute  in  equation  (69)  to  obtain: 

A3  - 146.457  kg  , T = C. 648566  x 1C~4  sec. 

c 

Continue  this  way  until  the  difference  between  two  subsequent 
3 

A or /and  two  subsequent  T is  sufficiently  small  to  be 

v 

negligible.  In  the  present  calculation  the  fourth  step  gives 
the  following  values: 

A3  - 146.763  kg 

Tc  = 0.645  x 10  4 sec. 

The  result  is  plotted  in  Fig.  18.  It  compares  very  well 
with  the  result  based  on  the  calculation  (c)  above. 

(ii)  Rectangular  ?lat4s  simply- supported  at  the  outer  edges : 


The  length,  a,  of  all  plates  in  x-direction  is  taken 
constant,  a = 20  cm.  The  calculations  are  then  made  for: 


a/b  =1  , a/b  =2  and  a/b  =4 . 


where  b is  the  length  of  the  other  edge.  Circular  frequency 
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a.„  are  given  by  equation  (36)  that  is: 
mn 

“on  ‘ “I0(m2  * “2  a2'b2)  ' “10  " *2/a2 

For  the  present  problem,  a^Q  ■ 0.622068  x 104sec.  Cal- 
culations are  carried  out  using  both  methods  with  the  results 
given  below: 


Square  plate  a/b  » 1 

F * sin  \ t (shape  a) 
Fx-  131.782  kg 

>»  ■ 4.636  x 104  1/sec 

Tc  * 0.678  x 10~4sec 


132.308  kg 
4.871  x 104 
0.645  x 10"4 


1/sec 


sec 


F - F1sin3^2  \t 
?1  - 146.763  kg 

\ * 4.875  1/sec 
Tc  - 0.644  x 10”4sec 


'I 


J 


(Generalized  Galerkin) 


(Collocation  (a)  ) 


(Collocation  (c)  ) 
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The  results  are  plotted  in  Fig.  19  next  to  the  results  of 
K.  Karas  [6],  based  on  the  reliable  step-by-step  integration 
method.  The  results  compare  very  well  with  those  of  K.  Karas. 
Shape  (c)  of  the  Collocation  method  can  not  be  distinguished 
from  the  Karas  result. 


Rectangular  plate  with  a/b  * 2 
F * F^sinVt 
?1  « 130.72  kg 
\ * 4.693  x 104  1/sec 

T = 0.669  x 10  4sec 

G 

Fx  * 132.680  kg 

\ • 4.829  x 104  1/sec 

Tq  * 0.650  x 10  4sec 

F * F1sin3^\t 
?1  - 145.967  kg 
\ • 4.844  x 104  1/sec 

Tc  - 0.649  x 10"4sec 


(Generalized  Galerkin) 

J 

^ (Collocation  (a)  ) 


» (Collocation  (c)  ) 


J 


See  also  Fig.  20 . 
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Rectanqular  plate  with  a/b  * 4. 

# » _ 

Only  the  functions  R^  , R£  , R^  # ^ an<*  ^ are 

plotted.  See  Figures  7,  and  14  to  17.  No  calculation  is  made 
for  the  contact  force. 

Circular  plate  clamped  at  the  outer  edge: 

Table  2 gives  a list  of  quantities  pma  and  cm  , the 
first  of  which  is  the  root  of  frequency  equation  (32)  and  the 
second  is  given  by  the  second  of  equations  (33) . The  procedure 
for  the  calculation  of  the  contact  force  is  identical  to  those 
given  above.  No  calculation  has  been  carried  out. 

Deflection 

The  maximum  deflection  is  computed  for  the  circular  plate 
and  square  plate  as  a function  of  time  and  plotted  in  Figs.  21 
and  22.  In  these  computations,  the  expressions  used  for  the 
contact  force  are  those  obtained  by  the  Generalized  Galerkin 
method,  for  the  circular  plate  and  collocation  (c)  for  the 
square  plate;  namely: 

F (t)  - 133.42735  sin  (4.627699  x 104t)  ♦ 1.505767 

sin  (9.255398  x 104t)  - 11.980813  sin  (13.883097  x 
104t)  (Circular  Plate) 

F (t)  - 146.763  sin3/2  4.875287  x 104t  * 134.334 
sin  (4.875287  x 104t)  - 14.926  sin  (14.625861 
x 104t)  kg  (Square  plate) 
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Deflection  obtained  by  use  of  a Dirac  8-function  ty$e  of  contact 

force,  having  the  same  impulse  as  the  contact  force,  placed  at 

the  mid-point  of  the  contact  interval  is  also  calculated.  That 

T 

is  F(t)  * 8 (t  - Tc/2)  with  F^  ■ ,£  F(t)dt  is  used  in  place 

of  F (t) . The  result  for  the  deflection  is  also  plotted  in  Fig.  22. 

Comparison  is  in  favor  of  all  the  above  shapes.  The  error 
committed  even  in  the  extreme  case,  where  the  8-function  is 
used,  seems  to  be  small. 
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Table  1 

pma'  cm  ' Simply  supported  circular  plates, (equation  29) 


m 

pa 

m 

m 

1 

2.2042 

1.7397 

2 

5.4451 

4.3309 

3 

8.6061 

6.7909 

4 

11.7578 

9.2557 

5 

14.9022 

11.7270 

6 

18.0438 

14.2004 

7 

21.2120 

16.6556 

8 

24.3536 

19.1249 

9 

27.4916 

21.5844 

10 

30.6332 

24.0505 

11 

33.7748 

26.5164 

12 

36.9164 

28.9909 

13 

40.0579 

31.4606 

14 

43.1995 

33.9271 

15 

46.3411 

36.3924 

16 

49.4827 

38.8631 

17 

52.6243 

41.3186 

18 

55.7659 

43.7099 

19 

58.9075 

46.2670 
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Table  2 

pffia,  cm  - Jlamped  circular  plates, (equation  32) 


m 

p_a 

m 

Cm 

1 

0 

0 

2 

3.1962 

2.7348 

3 

6.3064 

4.9073 

4 

9.4400 

7.3957 

5 

12.5771 

9.8695 

6 

15.7164 

12.3374 

7 

18.8565 

14.8025 

8 

21.9960 

17.3082 

9 

25.1376 

19.7519 

10 

28.2811 

22.1122 

11 

31.42C0 

24.6767 

12 

34.5613 

27.1385 

13 

37.7022 

29.7442 

14 

40.8438 

32.0801 

15 

43.9851 

34.5495 

16 

47.1266 

37.0096 

17 

50.2682 

39.4633 

18 

53.4097 

41.9252 

19 

56.5513 

44.3813 

20 

59.6929 

46.8340 
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Table  3 

Functions  R^,  Rg  for  the  Collocation  Method,  central 
impact  on  simply  supported  circular  plates,  (equations 58-60) 


Mo 

R1 

R2 

.040 

24.1985 

20.3397 

.080 

12.0880 

10.3370 

□ 121 

7.9639 

6.8764 

.160 

5.9109 

5.1983 

.200 

4.8442 

4.1920 

.240 

3.9143 

3.4747 

.300 

3.1871 

2.7443 

Functions 

Rl'  R2' 

Table  4 

Rp  R£  for  the  C 

ollocation  Method,  central 

impact  on 

simply  supported  square 

plates. 

a/b=lf  (equations 58 

Mo 

R1 

R2 

R1 

R2 

.08 

6.4559 

7.8811 

7.3765 

7.7704 

.12 

4.3901 

5.2337 

5.0066 

5.1442 

.15 

— 

— 

— 

3.9988 

.16 

3.3022 

3.8910 

3.7524 

3.8175 

.20 

2.6825 

3.2444 

3.0595 

3.2532 

.24 

2.2370 

2.5896 

2.5528 

2.5076 

.30 

1.7628 

1.9699 

1.9959 

1.8594 
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Table  5 

Functions  R^,  Rg,  RJ,  for  the  Collocation  Method,  central 
impact  on  simply  supported  rectangular  plates,  a/b*2, (equations 
58-60) 


**10 

R1 

R2 

R1 

R2 

.06 

4.3758 

5.2139 

5.0107 

5.0937 

.08 

3.2067 

3.7808 

3.6464 

3.6454 

.12 

2.1855 

2.7429 

2.4989 

2.8142 

.15 

1.8068 

2.3249 

2.0810 

2.4037 

.18 

1.5607 

2.0117 

1.8066 

2.0628 

o 

to 

o 

1.4267 

1.7990 

1 .6501 

1.8059 

.24 

1.2051 

1.4247 

1.3841 

1.3580 

. 3u 

0.9519 

1.0679 

1.0770 

0.9893 

Functions 

Rl'  R2* 

Table 

RJ.  R£  for  the  C 

6 

ollocation 

Method,  central 

impact  on 

siraoly  supported  rectangular  plates 

, a/b=4, (equations 

58-60) 

**10 

R1 

R2 

R1# 

V 

.06 

2.3161 

2.9355 

2.6751 

2.8738 

a 

o 

00 

1.7519 

1.9833 

1.9950 

1.8150 

.10 

1.3657 

1.3978 

1 . 5080 

1.1536 

.12 

1.0435 

0.8344 

1.1369 

0.5779 

.14 

0.8102 

0.5281 

0.8584 

0.2954 

.16 

0.6319 

0.2723 

0.6500 

0.0417 

.20 

0.3326 

0.0213 

0.3085 

-0.1410 

.24 

0.1620 

-0.0011 

0.2219 

-0.0590 

.32 

— 

0.0047 

— 

0.0141 
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Fig.  12  Functions  for  the  Collocation  Method 


Fig.  13  Functions  for  the  Collocation  Method, 

central  impact  on  simply  supported  circular  plate 
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i 


Fig„  14  Function*-  for  tie  Collocation  list  tied, 
central  impact  on  simply  supported  rectangular  plate,  (a/b=l,2,4) 


60 


i 


central  impact  on  simply  supported  rectangular  place,  (a/b-] 


125 


Fig.  18  G Gulerkin,  Collocation;  Forde-time  curves  for  central 
impact  on  a simply  supported  circular  plate 


150 


64 


l 


Fig.  19  G. Galerkin,  Collocation  (a),  (c) ; Force-time  curves 
for  central  impact  on  simply  supported  square  plate 


15C 


240 


Fig.  21  Deflection  curve,  circular  plate 
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